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Procedures for time-ordering the covariance function, as given in a previous paper (K. Kiyani and 
W.D. McComb Phys. Rev. E 70, 066303 (2004)), are extended and used to show that the response 
function associated at second order with the Kraichnan-Wyld perturbation series can be determined 
by a local (in wavenumber) energy balance. These time-ordering procedures also allow the two-time 
formulation to be reduced to time-independent form by means of exponential approximations and 
it is verified that the response equation does not have an infra-red divergence at infinite Reynolds 
number. Lastly, single-time Markovianised closure equations (stated in the previous paper |]J) are 
derived and shown to be compatible with the Kolmogorov distribution without the need to introduce 
an ad hoc constant. 
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I. INTRODUCTION 

In a previous paper 0, the Kraichnan-Wyld perturba- 
tion expansion was used to justify the introduction 
of a renormalized response function connecting two-point 
covariances at different times. The resulting relationship 
was specialized by suitable choice of initial conditions to 
the form of a fluctuation-dissipation relation (FDR). This 
was further developed to reconcile the time-symmetry of 
the covariance with the causality of the response func- 
tion by the introduction of time-ordering along with a 
counter-term. We pointed out that this formulation pro- 
vides a solution to an old problem in turbulence theory: 
that of representing the time-dependence of the covari- 
ance and response by exponential forms . We showed 
that the derivative (with respect to difference time) of 
the covariance with this time-ordering now vanishes at 
the origin. This allows one to study the relationships be- 
tween two-time spectral closures and time-independent 
theories such as the Fokker-Planck theory of Edwards 
or the more recent renormalization group approaches. 
We also showed that the renormalized response function 
is transitive with respect to intermediate times and re- 
ported a new Langevin-type equation for turbulence. 

In this paper we interpret the second-order response 
function as a mean-field propagator and show that in ad- 
dition to propagating two-time covariances it also links 
single-time covariances. We then make use of its newly 
established properties to re-derive the local-energy trans- 
fer (LET) response equation |tJ and show that it now 
contains a counter term which removes the singularity 
of previous propagator equations at t = t'. We also 
introduce a partial-propagator representation and hence 
reformulate the LET statistical equations. Furthermore 
we specialize the two-time equations to time-independent 
form by introducing exponential time dependences and 



show that the closure is well-behaved in the limit of infi- 
nite Reynolds number. Lastly, by Markovianizing time- 
history integrals, we end up with a Langevin-type theory 
which is compatible with the Kolmogorov spectrum with- 
out the need to introduce ad hoc constants as in the case 
of the EDQNM model for example. 

We begin by reviewing the subject of turbulence clo- 
sures and then go on to consider various aspects of ap- 
plying the FDR to noncquilibrium, macroscopic problems 
such as fluid turbulence. We begin by stating the basic 
equations. 



A. The basic equations 

Following standard practice in this topic ||}, we con- 
sider the solenoidal Navier-Stokes equation (NSE) in 
wavenumber (k) space, as follows: 

^ + v k 2 ^u a (k,t) =M a ^(k) J d 3 ju {i,t)u 7 {k-i,t) 

. (1) 

while the continuity equation for incompressible fluids is 

k a u a {k,t) = 0. (2) 

The inertial transfer operator M a p^(k) is given by 

M Q/37 (k) = {2i)- 1 [k p P ai {k) + fc 7 P Q/3 (k)], (3) 

while the projector P a p(k) is expressed in terms of the 
Kroneckcr delta as 



P a p{k) = 8 a p 



k a kp 

"ikP 



(4) 
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In order to introduce a statistical treatment, we shall 
denote the operation of performing an ensemble average 
by angle brackets, thus (• • •), and restrict our attention 
to isotropic, homogeneous turbulence, with energy dissi- 
pation rate e and zero mean velocity. As a result of this 
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restriction, the covariance of the fluctuating velocity field 
takes the form 

(u a (k, t) Uf3 {k', t')) = C{k- 1, t')P aP (k)S(k + k'), (5) 

where a, ft = 1,2 or 3. The corresponding single-time 
quantity may be written as 



C(k;t,t) = C(k,t), 



(6) 



where the single-time one-point covariance C(k,t) may 
be interpreted as a spectral density and is related to the 
energy spectrum by 



E(k,t) = 4Trk 2 C(k,t). 



(7) 



Using l(T|l. we can also obtain an equation describing the 
energy balance between spatial modes. To do this, we 
first multiply each term in Q by u a (— k,t). Then we 
form a second equation from Q for it ff (— k, t), multiply 
this by u Q (k, i), add the two resulting equations together, 
integrate over k' and average the final expression. This 
leaves us with 



d 



f t +2v k 2 ) P QCT (k)C(fc,i) 



M al3l (k) J d 3 jC <J/ 5 7 (-k,j,k - j;t) 
-M CT/37 (k) y d 3 jC CT/37 (k,j,-k-j;/), (8) 



where 



Ca/3 7 (k,j, -k- j;t) = (ii Q (k,t)w /3 (j,i)u 7 (-k - j,t)), 

(9) 

and where we have also used the property 

M a/37 (-k) = -M a/37 (k). (10) 

If we then take the trace of by setting a = a and 
summing over a (noting that Tr P a p = 2) , and multiply 
each term in (jSJ by 27rfc 2 , we obtain 



li + v k 2 j El/.-./) 7( /,-./). (I!) 



where 

T(M) = 27rfc 2 M Q/37 (k) y"d 3 j {C Q/37 (-k,jk-j,/) 

-C Q/37 (k,j,-k-j,/)}. (12) 

Evidently, in order to solve for the energy spectrum (or, 
second-order moment) we need to know the third-order 
moment. Hence we are faced with a hierarchy of statis- 
tical equations to be solved; and this is, of course, the 
notorious closure problem. 



B. Eulerian statistical closures for isotropic 
turbulence 



In order to study isotropic turbulence, we have to add 
a noise term or stirring force to the right hand side of the 
NSE, as given by QJ. Denoting this term by / Q (k, t), we 
specify it in terms of its distribution, which we take to 
be Gaussian, and its covariance, which we take to be of 
the form 

</„(k, t)f (k', t 1 )) = W(k)(2ir) 3 P aP (k)6(k + k')S(t + t'). 

(13) 

We note that W(k) is a measure of the rate at which the 
stirring forces do work on the fluid and for stationarity 
must satisfy the condition 



Ank 2 W(k)dk 



2vk 2 E(k)dk. (14) 



The perturbative treatment of the equations of motion is 
based on an expansion about a Gaussian zero-order veloc- 
ity obtained by solving the NSE with the nonlinear term 
set to zero. The resulting expansion shows clearly 0] the 
effect of nonlinear mixing such that any correction to the 
zero-order field must have a non-Gaussian distribution, 
which indeed is implied by the existence of the third- 
order moment (and the existence of inter-modal energy 
transfer) . Renormalization of the perturbation expansion 
corresponds to either partial summation or term-by-term 
reversion: for details reference should be made to the pa- 
per by Wyld 3] and the books by McComb |9j and Leslie 
|5j . Our present interest is restricted to the second-order 
equation for the velocity covariance, which is obtained by 
this procedure, thus: 



d_ 



vk 1 



C{k;t,t') 



J d 3 jL(k,j) 



dsR(k;t', S )C(j;t, S )C(\k-j\;t, S ) 



dsR(j;t,s)C{k;s,t')C(\k-]\;t,s) , (15) 
and on the time diagonal 

(Jt + 2vk 2 ^ C(k,t) = 2 J d 3 jL (k, j) 

x / dsR(k;t, s)R(j;t, s)R(\k — j| ;t, s) 
Jo 

x [C (j, s) C (|k - j| , s) - C (k, s) C (|k - j| , s)} , (16) 

where the coefficient L(k,j) is given by: 

L(k,j) = -2M a ^{k)M 0a5 {i)P l5 {k- j). (17) 

This may be evaluated in terms of the scalar magnitudes 
k, j and /i = cos 9, where 9 is the angle between the two 
wavevectors k and j, thus: 

L(k j) - [Mfc2 + ]2) kj(1 + V)] ^ 1)kj (18) 
k 2 + j 2 — 2kjfi 
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It should also be noted that the coefficient L(k, j) is sym- 
metric under interchange of the two wavevectors: we shall 
use this fact presently to demonstrate conservation of en- 
ergy. 

At this stage we should note that for this to be a 
closed set of equations for the covariance C, one has to 
have an additional equation to determine the response 
function R. Equation l|15fl wa s originally derived by 
Kraichan. This closure was completed by an equation for 
the response-function R(k; t, t'), and is known as the di- 
rect interaction approximation or DIA. The basic ansatz 
of DIA is that there exists a response function such that 



<5u a (k, t) 



R a0 (k;t,t')6f (k,t')dt', (19) 



and that this infinitesimal response function can be 
renormalized. The resulting response equation is 



d_ 

dt 



R(k;t,t') + J d 3 jL(k,j) 



x J dt"R(k;t",t')R(j;t,t")C(\k-j\;t,t") 

= 6(t-t'). (20) 

Later Edwards derived a time-independent covariance 
equation by the self-consistent introduction of a gener- 
alized Fokker-Planck equation as an approximation to 
the (rigorous) Liouville equation. We shall refer to this 
theory as EFP, and this along with the more general self- 
consistent field (SCF) theory of Herring and the DIA 
make up our trio of pioneering spectral closures. Further 
discussion can be found in the books ^ n the lit- 

erature, much attention has been given to the fact that, 
although these theories have many satisfactory features, 
they are all incompatible with the Kolmogorov (K41) 
power law for the energy spectrum E(k) pT| . However, 
in the present paper we shall concentrate on only a few 
key points. The first of these is that the covariance equa- 
tion of the EFP theory can be shown to be equivalent, for 
the stationary case, to the second-order truncation of the 
Kraichnan- Wyld perturbation theory, if we assume expo- 
nential time dependences. That is, the EFP covariance 
equation can be obtained by substituting into equation 
(|15fl the following assumed time dependences: 



C(k, t-t') = C(k) exp{-uj(k)\t - t'\}; (21) 



and 



R(k,t-t') = exp{-u>(k)(t - 1')} for t > t'; 

= for t< t'. (22) 

Then, integrating the right hand side of equation (I15|) 
over intermediate times, one obtains (with some re- 
arrangement): 



W(k)-2vk 2 C(k) = / d 3 jL(k,j) 



C(|k-j|)[C(fc)-C(j)] 
w(*)+w(j)+w(k-j|)' 
(23) 



where we have added the term W(k) to the energy bal- 
ance in order to sustain the turbulence against viscous 
dissipation. Equation (|23|l is just the form originally de- 
rived by Edwards [g- 

This simple form is helpful in understanding certain 
properties, such as the conservation of energy by the non- 
linear term and the behaviour of the system in the limit 
of infinite Reynolds number. For instance, integrating 
both sides of equation (|23|l with respect to k and invok- 
ing equation ljTI|) leads to: 



0. 



(24) 



where the vanishing of the right hand side results from 
the antisymmetry of the integrand under interchange of k 
and j. This result helps us to interpret the EFP response 
function or eddy decay rate R(k) which takes the form 



R(k) 



d 3 jL(kJ] 



C(|k-j| 



w(fc)+w(i)+w(k-j|)' 



(25) 



At the time this was interpreted as allowing one to write 
the energy balance equation as: 



W(k) - 2vk 2 C(k) = R{k)C(k) 

,3,r^ * C(|k-j|)C(j) 



d 6 jL{k,j) 



w(k)+w(j)+wQt-j\y 



(26) 



or, the eddy decay rate represented the loss of energy 
from mode k due to energy transfers to all other modes. 
The situation is more complicated for DIA, but the anal- 
ogous comment has been made by Kraichnan Q that the 
energy loss from mode fc is directly proportional to the 
excitation of that mode, viz., C{k]t,t'). 

Later, it was pointed out that an ad hoc modification 
could be made to EFP by noting that the entire energy 
transfer term (i.e. the right hand side of acts as an 
energy loss in some regions of wavenumber space whereas 
in others it behaves as an input. This led to a definition of 
the response which was compatible with the Kolmogorov 
spectrum [l2l Il3| and this was subsequently generalised 
to the two-time local energy transfer or LET theory . 



C. Fluctuation-dissipation relations (FDR) 

It is well known that the response of microscopic sys- 
tems in thermal equilibrium to small perturbations is 
fully determined by the covariance of the system fluc- 
tuations about equilibrium. In our present notation, the 
relationship may be written as 



C(k;t,t') = R(k;t't')C(k;t',t'), 



(27) 



which is the fluctuation- dissipation relation or FDR. This 
result was extended by Kraichnan to nonlinear dynamical 
systems in thermal equilibrium |l5l Il6j | and by Leith [T7| 
to inviscid two-dimensional chaotic flow. Also, Deker and 
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Haake ,18j give several examples of classical processes for 
which a FDR will hold and these include (of particular 
relevance to the present discussion) forced viscous flows 
where the stationary probability distribution is Gaussian. 
In realistic cases, such flows will have a non-Gaussian 
distribution due to nonlinear mode coupling. However, 
one case of interest arises in a pioneering application of 
renormalization group methods to stirred fluid motion 
[l9| , where a fluctuation-dissipation relationship is found 
to hold in the limit k — > 0. That is, the long-wavelength 
behaviour at lowest nontrivial order of perturbation the- 
ory. 

The paper by Leith is particularly interesting. While 
recognising that the FDR cannot apply exactly to real 
fluid turbulence, it puts forward rather convincing heuris- 
tic arguments for believing that it could be a reasonable 
approximation. Also, it cites the investigation of Her- 
ring and Kraichnan [20| in support of this view. Here 
the non-stationary generalization of the SCF ,21], which 
differs only from DIA in the use of the FDR, gives very 
similar results to it. We shall discuss this use of FDR in 
more detail later, when we consider its role in the LET 
theory. 

Lcith's optimistic view not only inspired successful 
practica l ap plications of the FDR to study climate sen- 
sitivity |22j ] and viscosity renormalization pflj . but was 
also seen as seminal in stimulating an important series 
of papers which examined the applicability of the FDR 
from the point of view of dynamical systems theory |24j - 
|26j . The overall conclusion of these papers can be stated 
as follows: 

1. A general relationship exists for the response of a 
chaotic system in terms of its stationary probabil- 
ity distribution provided that the system is dynam- 
ically mixing. 

2. If the stationary probability is Gaussian in form, 
then the relationship reduces to the FDR as given 
by equation (|27|l . 

Of course in real fluid turbulence the probability distribu- 
tion is not Gaussian, nor is it known exactly. However, 
as we have shown in 0, the FDR can be derived for 
turbulence to second-order in renormalized perturbation 
theory and hence, if used appropriately, is consistent with 
a closure approximation of that order. We shall return 
to this point later. 



D. The time-ordered FDR 

In pj we postulated that in the context of the 
Kraichnan- Wyld perturbation theory we may re-write 
the existing relationship between the zero-order covari- 
ance and zero-order response in a renormalized form as 

C aa (k; t, = (t - s) i? Q£ (k; t, s) C ta (k; a, t') , (28) 



or in its isotropic version as 

C (k; t, t') =9{t-s)R [k\ t, s) C (k; s, t') , (29) 

where the Heaviside unit-step function 9(t — s) explicitly 
states the causality condition. As yet we have taken no 
decision about the ordering of the two times t and t', 
and thus the symmetry under interchange of t and t' is 
untested in 

If we explicitly state the time ordering as t > t' say, 
then this is equivalent to applying 9(t — t') to both sides 
of 

9(t- t') C (k; t, t') =9(t-t')9(t-s)R (fc; t, s) C (k; s, t') , 

, (30) 

and this is the beginning of the LET theory. In it, we 
have postulated the existence of a renormalized propaga- 
tor and have made use of the Heaviside unit-step function 
to make the time-ordering manifest. 

The generalized fluctuation dissipation relationship is 
obtained by setting s = t' in (|3TJ1) to get 

9(t- t') C (k; t, t') =9{t-t')R {k; t, t') C (k; t', t') , 

where the time-ordering is set by the requirement s = t' . 

In Q we introduced a representation of the covariance 
which preserves its symmetry under interchange of time 
arguments, thus: 

C(k;t,t') = 9{t-t')C{k;t,t') + 9{t' -t)C{k\t,t') 
-8 t ,t>C (k;t,t') . (32) 

We can easily show that this representation does what 
it is supposed to do by looking in turn at the separate 
cases: t < t'; t > t' and t — t' , and this is left for the 
reader. 

Now, using H30fl to expand the right hand side of l|32() 
we obtain 

C(k;t,t') = 9(t-t')9(t-s)R(k;t,s)C(k;s,t') 

+9 (t' -t)9 (t' -p)R(k;t',p)C(k;p, t) 
-St, v C{k; t,t'). (33) 

Or, this result may be written more like the FDR by 
instead using l|31|l to construct it 

C(k;t,t') = 9(t-t')R(k;t,t')C{k;t',t') 
+9 (t' - t) R(k;t' ,t) C (k;t,t) 
-6t,fC (k;t,t') . (34) 

The symmetry of both these covariances, (|33f) and 
can be broken by applying a unit-step function to both 
sides. This will yield either l|3(J|) or (|31|l . depending on 
which time-ordering we choose. 

II. THE PROPERTIES OF THE MEAN-FIELD 
PROPAGATOR 

In this section we begin by reviewing the introduction 
of a velocity propagator, as in the original formulation 
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of the LET theory [3 and note that in this context the 
propagator introduced in Q] is a mean-field propagator. 



A. The velocity field propagator 

From the exact solution of the solenoidal NSE (see 
we have: 



(k,t) = R^(k;t,s)u a (k,s) + 
+ \\ /Vi?(°J(k;M'') 



x / rf 3 jM^ 7 (k)^(j,f> 7 (k-j,t") 



(35) 



where is the 'viscous' or zero-order response tensor 
and the hat is used to emphasize that this is the 're- 
sponse' associated with the instantaneous velocity field. 

Expanding u a (k, t) in a perturbation series around a 
gaussian solution and equating zero-order terms we can 
say that the equality 



«i 0) (M)=i^°j(k;MK 0) (k,. 



(36) 



illustrates the propagator-like nature of Raa (k;£, s). 
Then from looking at the form of l(35() . we can postulate 
the existence of a renormalized propagator such that we 
obtain a renormalized version of ((36(1 



u Q (k, t) = R aa (k; t, s) u a (k, s) . 
Multiply 1(23 by tt« (-k,f) 



(37) 



u a (k, i) up (-k, i') = i? Q(T (k; t, s) w CT (k, s) (-k, i') , 

(38) 

and average this equation to obtain 



C a p(k; t, t') = R aa (k; t, s) C a p (k; s, t') , 



(39) 



where the propagator is statistically independent of the 
velocity field and we have used the mean-field approxi- 
mation 

(i? aCT (k;t,s)J> = R aa (k;t,s). (40) 

As usual, equation can be turned into a simpler 
scalar form by using the properties of isotropic tensors 

C (k; t, t') = 9{t-s)R (k; t, s) C (fc; s, t') . (41) 

The transitivity of R ar (k;t, s) with respect to inter- 
mediate time can be proved by applying equation (|3Tjl to 
the right-hand side of itself 

u a (k, t) = R aa (k; t, s) R ap (k; s, t') u p (k, t') , (42) 

and realising that we could also have written this as 

u a (k, t) = R ap (k; t, t') u p (k, £') , (43) 

implying the result: 

R ap (k; t, t') = R aa (k; t, s) R ap (k; s, t') , (44) 

and t > s > t'. 



B. The mean- field propagator 

The simple property of the propagator 
R(k;t,t) = 1, 



(45) 



can be easily shown to be necessary by setting s = t in 
(l3THl . Other properties can be obtained by equating the 
right hand side of H3()|l with the right hand side of 131(1 : 

0{t-t')R(k;t,t')C(k;t',t') 

= 6(t-t')6(t- s)R (fc; t, s) C (fc; s, t') . (46) 

Expanding the right hand side of ||S^ using we ob- 
tain 



9(t-t')R(k;t,t')C(k;t',t') 

= [9(t-t')9(t- s)R(k;t,s)x 
x6(s-t')R(k;s,t')C(k;t',t')} > " 

+ [6{t-t')6 (t- s)R(k;t,s) x k , 
xB(t'- s)R{k;t',s)C{k;s 7 s)} > " 

-[6{t-t')d(t-s)R(k;t,s) ' 
x S t , )S C(k;s,t')}. 



(47) 



Dividing the right hand side into three groups of terms 
labelled respectively a, b and c, we will now look at 147(1 
for the two separate cases:- Case 1, t > s > t'; and 
Case 2 t > t' > s. 



1. Transitivity with respect to intermediate times 

Here we consider case 1 corresponding to t > s > t' 
and implying that 6 = and c = in equation ((47(1 . This 
will leave 

6(t-t')R(k;t,t')C(k;t',t') 
= [6(t-t')9(t- s)R(k;t,s) x 
x9(s-t')R (k; s, t') C (k; t' , t 1 )] . (48) 

We now use the contraction property of the Heaviside 
function 

9(t- s) 9 (s -t')=9(t- t') , (49) 

to write |@S) as 

9(t- t') R(k;t,t') C (k; t', t') 

= 9(t-t') R (fc; t, s) R (fc; s, t') C (fc; t' , t') . (50) 

From this above result, we can deduce the transitive 
property of the propagator 



R{k;t,t') =R(k;t,s)R(k;s,t') 



(51) 



This result also tells us that the transitivity of the prop- 
agator holds only for times s which are intermediate be- 
tween the two times t and t! . This makes sense because 
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otherwise, if s was outside the range between t and t', 
we would have propagation backwards in time which vi- 
olates causalit y. This is a result which was previously 
only assumed [27|| on the basis that it could be ex- 
pected to follow from the corresponding relationship for 
the velocity-field propagator, and is now proved. 



2. Linked single-time covariances 

Next we consider Case 2 t > t' > s, which corre- 
sponds to a — and c = 0, leaving 



6{t-t')R(k;t, t') C(k;t',t') 
= [8 {t-t')6(t- s)R(k;t,s) 
x 6{t' - s)R(k;t',s)C(k;s,s) 



(52) 



This result is important because it links two single-time 
covariances. This fact becomes clearer if we take the 
special case of t = t'. This gives 

C (k; t,t) = 6 (t — s) R (k; t, s) R (k; t, s) C (k; s, s) , 

(53) 

implying that we need two propagators to link single-time 
covariances. Defining 



R(k;t,s) :=R(k;t,s)R(k;t,s), 



(54) 



Equation (|53|l can be modified to make it look like equa- 
tion Q 



C (k, t) = 6(t-s)R (k; t, s) C (k, s) 



(55) 



Again, the presence of the unit-step function ensures that 
the covariance can only propagate forwards in time. 



III. DERIVATION OF THE LOCAL ENERGY 
TRANSFER (LET) RESPONSE EQUATION 

The starting point for the LET theory is the second- 
order renormalized covariance equation as given by (|15|l . 
We can now proceed in two ways. 

1. The first is to substitute (|34fl in l|15|) and then 
choose t > t'. 

2. The second is to choose t > t' and multiply both 
sides of (I15fl by 8 (t — t') to show the range over 
which the equation will be valid. Then follow this 
by using the FDR, in the form of equation (|31fl . 
throughout. 

Note it is important that we do not set t = t 
in the covariance equation (|15|) as we can only 
do this after evaluating the derivative of the 
2-time covariance. 



Both methods are equivalent but the second is the easier 
to use in practice. Thus we begin by choosing the time- 
ordering to be t > t' and multiplying (|15|l by 0(t — t') 

6(t - t')^-C (k; t, t') + 9(t - t')vk 2 C (k; t, t') 
at 



6{t-t') J d 3 jL(k,j) 



x <^ / dsR(k;t',s)C(j;t, S )C(\k-j\;t,s) 



f dsR (j; t, s) C (k; s, t') C (|k - j| ; t, s)\ . (56) 



Let us look at the first term of the left hand side of 
e(t-t')-^C(k;t,t') 



nt J (t - t') C (k; t, f) - C (k; t, f) ^8 (t - f) 



d_ 

d 

m 1 



(t-t')R(k;t,t')C{k;t',t') 



-C(k;t,t')^- t 0(t-t') 



(57) 



where we have applied the product rule in the second line, 
and the FDT, lpTT|) in the third line. After substituting 
the differential of the Heaviside unit-step function 



8 

0t< 



■{t-t') = 6(t-t'), 



(58) 



we reach our final form for this part of the response equa- 
tion. Thus: 

left hand side of (|56|l 

'' "(t-t')R(k;t,t')C(k;t',t')-C(k;t,t')5(t-t') 



at 

+vk 2 9 (t - t') R (k; t, t') C (k; t' , t') , 



(59) 



where the FDR was used on the second term of the 
left hand side of (|56|l also. 

Now we evaluate the second time integral on the right 
hand side of which we label as TI 2 : 

Th = 6{t-t') f dsR(j;t,s)C(k;s,t')C(\k-j\;t,s). 
Jo 

(60) 

We need to have the appropriate 8 functions in front of 
the covariance so that the broken time-reversal symmetry 
becomes manifest. This information is present in the 
arguments of the propagator and in 9(t — t'). So for 
C(|k-j|;M) 



8{t-t') f dsC(\h-j\;t,s) 
Jo 



= 8 {t-t') / ds6(t-s)C(\k-j\;t,s), (61) 
Jo 
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and for C (k; s, t') 



(t-t') / dsC(k;s,t') 



= 6{t-t') f dsC(k;s,t') 
Jo 

+6 (t-t') J dsC(k;s,t') 

= 6 (t-t') f ds9(t' -s)C(k;t',s) 
Jo 

+6 (t - t') J ds6 (s - t') C (k; s, t') , 



(62) 



where we have used the property C(k;t,t') — C(k;t',t) 
in the fourth line. With these results we can now write 

as 



Tin 



6 (t-t') dsR(j;t,s)9(s-t')x 



xC(k;s,t')6(t- s)C(\k-i\;t,s) 



(t-t 1 ) / dsR(j;t,s)6(t' - s) x 



x C(k;t',s)9(t- s)C(\k-j\;t,s) 



(63) 



The evaluation of the first integral on the right hand 
side of H56fl follows similarly so that the final LET re- 
sponse equation is given by 



(t - t') R (k- 1, t') C (k; t', t') - C (fc; t, t') S (t - t') 



+vk 2 9 (t - t') R (k; t, t') C (k; t' , t') 
= J d 3 jL(k,j)e(t-t') x 

dsR(k;t',s)9(t - s) x 



xC(j;t,s)9(t-s)C(\k-j\;t,s) 
t' 

dsR (j;t,s)6(t' - s) x 



" 



x C(k;t',s)6(t- s)C(\k-i\;t,s) 



dsR (j;t,s)6(s-t r ) x 



Jv 



xC(k;s,t')9(t-s)C(\k-j\;t,s)\} 



(64) 



Multiplying both sides by 9(t-t'), dividing by C (k; t' : t') 
and noting that 



6(t - f)C (k;t,f) 
C(k;t',t') 



9(t-t')R(k;t,t'), (65) 



from the FDR, in the form of equation (|31|l , we reach the 
simplified form with the broken time-reversal symmetry 
manifest 

9 (t - t') (J^ + vk 2 ^j 9(t- t') R (k; t, t') 



-9(t-t')R(k;t,t')5(t-t') + 



J d 3 jL(k 7 j)0(t-t') 



x / dsR(j;t,s)R(k;s,t')9(t- s)C (\k-j\;t,s) 



d 3 jL(k,j)6(t-t') / ds 



■(t-s)C(\k-j\;t,s) 



C(k;t',t') 



x {R(k;t',s)9(t - s)C(j;t,s) 
- R(j;t,s)9(t' -s)C(k;t',s)}. 



(66) 



A. Comparison with previous forms 

Apart from the addition of the second term on the left 
hand side 



9(t~t l )R(k- : t,t')5(t-t l ) , 



(67) 



equation l|66|l is the same as the LET response equation 
which appears as equation (3.19) in |7|, equation (20) in 
[2l| and equation (7.146) in 0- The natural addition of 
this extra term as a consequence of time-ordering, fixes 
the problem of the singularity in the time-derivative of 
the response equation (|66|l which occurs when one takes 
t = t'. More important, if we compare (64) with the DIA 
response equation (20), the additional terms on the right 
hand side of (64) cancel the infra-red divergence and en- 
sure compatibility with the Kolmogorov K41 spectrum. 



IV. THE TWO-TIME LET THEORY 

A. Partial-propagator representation 

We may write the propagator in a representation which 
separates the discontinuous part as a Heaviside unit step 
function, thus: 



R(k;t,t') = 9(t-t')TZ(k;t,t'), 



(68) 



where lZ(k; t, t') is a representation of the propagator but 
without the discontinuity at t = t'. So using l|6^|) and the 
FDR 127(1 to turn two-time covariances into single-time 
form, (64) for the response function becomes: 



6(t-t') (^- + iyk 2 ^TZ(k-t,t') 

= -9 (t-t') J d 3 jL (k,j) J dsTZ(k;s,t') 
xK(j;t,s)n(\k-i\;t, s) C (\k - j| , s) 
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+0 (t - t') j d 3 jL (k, j) j ds[n (k; f , s) n (J; t, s) 



xTZ (|k - j| ; t, s) C{ c {k ^ ) S) \C (j, s)-C (k, s)] } , (69) 

for t > t' . The counter-term has been cancelled by use 
of the product rule in the time-derivative. 



B. The LET closure equations 

The LET Equations may now be summarised as fol- 
lows. For the two-time covariance, we have equation 1)150. 



j t • i'lr ) C (/, :/./' 



J d 3 jL(k,j) x 
dsR (k; t', s) C (j; t, s) C (|k - j| ; t, s) 



dsR (j; t, s) C (k; s, t') C (|k - j| ; i, s) j , (70) 

and likewise equation (|16|) for the single-time covariance: 

(jt + 2vk^ C(k,t) = 2j d 3 jL (k, j) x 

x / dsR (k; t, s) R (j; t, s) R (|k - j| ;t,s) X 
Jo 

[C (j, s) C(|k - j| , s) - C (k, s) C (|k - j| , s)} . (71) 

where we have invoked the FDR so that all two-time 
covariances are turned into one-time covariances. 

For the response function we can use either (64) or 
(67). The above equations along with the generalised 
fluctuation-dissipation relation (FDR), 

6{t - t')C(k; t, t') = 0(t - t')R(k; t, t')C{k, t'), (72) 

from which the LET is derived, and the single-time co- 
variance link equation 



C(/c, t) = 0(t - s)R(k; t, s)R{k; t, s)C(k, s), 



(73) 



complete the set of LET equations. 

The LET equations have been applied, along with 
those of the DIA, to the problem of free decay of isotropic 
turbulence from arbitrary initial conditions, over a wide 
range of Taylor- Reynolds numbers > • I n these in- 
vestigations, the covariance equations on and off the time 
diagonal were solved simultaneously with the relevant re- 
sponse equation. It was later realized that for the LET 
theory, the response equation could be replaced by the 
FDR, as given by (70), and this reduced the computa- 
tional effort well below that of DIA: see 0, ISJ. This 
work was for three-dimensional turbulence, while an ex- 
tensive investigation of the two-dimensional case has also 
been carried out for DIA, SCF and LET theories |3lj . 

m, M. 



C. Behaviour in the limit of infinite Reynolds 
numbers 



The later two-time versions of the LET theory claimed 
that their solutions were compatible with K41. However, 
this was never shown explicitly. Compatibility with K41 
is now demonstrated for the LET response/propagator 
equation as given by 1)69(1 . We begin by writing (|69|l in 
stationary form. This means that all (single-time) co- 
variances become time independent: 

C(M)->C(fc), (74) 

and we write the propagator in relative time coordinates 

Tl(k;t,t') =TZ(k;t-t'). (75) 

Next we assume the exponential form for the propagator 

K(k;t-t') = exp[-u(k)(t - t% (76) 

where, as before, u>(k) is the total eddy-decay rate. These 
changes result in the response equation becoming: 



>(t-t 



') ^ + ^cxp[-c(fc)(t-0] 

-0(t-O J d 3 jL (k,j) C (|k - j|) J dsx 

x exp[-w(fc)(fl - - w(i)(* - «) - w(|k - - s)} 

t' 



6(t-t' 


J d 3 jL(k,j)^ 







x { exp[-w(k){t' -s)- w(j)(t -s)- w(|k - - a)] 

(77) 



x C( L k ~ j|) [C(j)-C(k)] 



C{k) 

Doing the differentiation, setting t = t' and carrying out 
the time integration results in an equation for u(k): 



u>(k) = vk 2 



d 3 ^(k,j) 



C(\k-j\)[C(j)-C(k)] 



C(k)[u(k)+w(j) + u(\k-j\)] 
x {l-exp[-(u(k)+u(j)+u(\^-3\))t])}, (78) 

where one can ignore the last term involving the expo- 
nential as we are considering stationary systems which 
are time independent. Another way to justify the ne- 
glect of this term is to realise that it originates from the 
fact that we chose to have the initial conditions at t = 
rather than the more usual t = — oo. 

To show that (|78|l is not divergent we complete 
our analysis by substituting the inertial range/ infinite 
Reynolds number forms for C{k) and w(fc): 



C{k) 



a£2/3 /; -H/3 



(79) 
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(80) the response equation 



where a and j3 are constants, and by writing the integral 
in k,j,fi variables 



(->U I') [ ■ vk? ) 7C (A : / 



cj(k) = vk 2 



, dli kj 3 (n 2 - l)Hk 2 +j 2 ) kj(l + 2 M 2 )] y 



(81) 



0(t-t') / d 3 jL (k,j) 



fc 2 + j 2 - 2kj/j, 

p-^/ 3 ik - jr ii/3 [j- u / 3 - k~ ii / 3 \ 



k -ll/3 



fc 2 /3+j2/3 + | k _j 



,2/3 



where /i is the cosine of the angle between the two vectors 
k and j. 

There are three possible sources of divergence (of the 
infra-red type) in this expression. However, from equa- 
tion H81|l. it may be seen that the k — > and j ' — > 
limits do not pose a problem. The final possible source 
of trouble |k — j| — > can be resolved by realising that 
the term [j~ n / 3 — fc~ n / 3 ] cancels the divergence caused 



by the |k-j| 



-11/3 



term. This is shown by expanding 



-11/3 



(k 2 +j 2 -2kjfi) 



-11/6 



(82) 



and substituting in equation 1)81 [I. One then Taylor ex- 
pands k around j to leading order in e = fc — j in both the 
numerator and denominator of the integrand in equation 
(|81|l . This results in the integrand becoming 



(ll/3)a/3 



-V/3 



(2)17/6 

y [2j 2 mi-m)-j 2 ] Qi 2 -i)(i-M)- 17 / 6 6 

X jW/6 [ 2j 2/3 + (2j2)l/3(l _ M )l/3] 



(83) 



and focusing on the term (/i 2 — 1)(1 — /i)~ 17 / 6 e we can 
see that as e — > 0, the integrand goes to zero, except at 
fx = 1 where the integrand is singular. This singularity 
can be avoided if we write the limits of the /i integral as 



1i 



(84) 



where f 1 implies in the limit approaching 1 from below. 

This completes the analysis in the limit of infinite 
Reynolds number. Further information on the above 
technique can be found in |(| . 



V. SINGLE-TIME MARKOVIANIZED LET 
THEORY 

The relevant single time LET equations are the single- 
time covariance (using the partial propagator form): 

(j t + 2vk 2 ^j C(k,t) = 2j d 3 jL (k, j) x 

x / ds1l(k;t,s)n(j;t,s)1l(\k-j\;t ) s)x 
Jo 

[C(j, S )C(|k-j|, s )-C(fc, S )C(|k-j| )S )], (85) 



dsTZ{k;s,t') x 



xK{j;t,s)K(\k-j\;t,s)C(\k-}\,s) 
+9(t-t') Jd 3 jL(k,j)J^ ds{n(k;t',s)1l(j;t,s)x 
x^(|k- j| ; t, s) C ^~i: s) [C (j, s)-C (fc, s)}\ , (86) 



C(k,f) 

and the single-time covariance link equation 

C(fc, t) = 9{t - s)R(k; t, s)R(k; t, s)C(k, s), (87) 
Making a Markovian approximation we can write i|85|) 



as: 



^+2!/fc 2 )c(M) = 2 J d 3 jL(k,j)D(k,j,\k-j\;t) 
xC(\k-j\,t)[C(j,t)-C(k,t)}, (88) 



where the Markovian approximation amounts to updat- 
ing each C(s) to C{t), and where 

D(k,j,\k-i\;t) 



dsTl{k;t,s)TZ(j]t,s)n(\k-j\;t,s), (89) 



is the memory time. 

We now need some way of computing -D(fc, j, |k — j | ; t) ; 
that is, of updating it. We do this by differentiating l(8H 
with respect to t 



d 



£>(fc,i,|k-j|;t) 



ds 



at 
= i 



+K(k;t,s) ( — K(j;t,s) ) ft (|k - j| ; t, s) 



-n{k;t,s)n{j;t,s) [ — n(\k-j\;t,s) 



^K(k;t,s) j ■Rij:1.s)R[ k-j :/.,) 



(90) 



To evaluate (|90|l we need to know the dynamical be- 
haviour of 1Z (fc; t, s). We obtain this from equation 18011 . 
Proceed by writing 180|) in Langevin form 



6{t - t') 



+ vk 2 + rj(k;t,t r ) 



Tl(k;t,t') = 0, (91) 



where 
V(k;t,t') 

= 0{t-t') J d 3 jL (k,j) J dsTZ(k;s,t') 
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xK(j;t,s)K(\k-j\;t,s)C(\k-j\,s) 
~6 (t ~ t') J d 3 jL (k, j) J ds[n (k; if, s) K (j; t, s) x 
xllQk - j| ; t, s) C{ c { ~^ S) [C (j, s)-C (k, S )]| , (92) 

is the turbulent eddy-decay rate and is obtained by com- 
parison with |(HEll • 

Rearranging l|91|) we obtain 

9{t-t')^- t n{k;t,t') 

= -6{t - if) [vk 2 + rj{k- 1, if)] K (k; t, if) , (93) 
and this allows us to write l|9()|l as: 



-D(k,j,\k-i\;t) 



ds 11 (k; t, s) K (j; t, s) K (|k - j| ; t, s) x 
(vk 2 + vj 2 + v \k - j\ 2 ^j + rj(k;t,s) 



+ ?7(j;M) +»7(|k-j|;i,s) 



(94) 



To be able to calculate (1941 we need to take the Marko- 
vian step 



r)(k;t,s) -> r)(k,t). 



(95) 



We can justify this step by looking at equations (|57|) . 
and l|91|l . Equation (|91|l has a general solution 



K (k; t, t') = exp <^ -vk 2 (t -t')- \ ds7j(k; s, t') 



(96) 

If we write (I88fl in the suggestive form 

C(k,t) = 0, (97) 



^ + 2^fc 2 + 2e(fc,t) 
ot 



where 



i(k, t) = - j cfjL (k,j) D(k,j, |k - j| ;t) X 

then we can write the general solution of (|97fl as 

C(k,t) 

= exp|-2i//c 2 (t-t') -2 J ds£(k,s)Xc(k,t') 

exp j-i^ 2 (i -t')- J t ds Z( k , s )j C (*> ■ (99) 
Writing as 

C(k, t) = 6(t - t')Tl{k; t, t')1Z(k; t, t')C(k, t'), (100) 



and comparing with (|99(l , this suggests 
K (fc; t, t') = exp l^-vk 2 (t - t') - jT cfe£(fc, s) j . (101) 
But comparing this with l|96|) . we see that 

ds£(k,s) = / dsr](k;s,t r ). (102) 



Comparing the forms of £(fc, s) and rj(k;t,s), (|98fl and 
(IHll, we see that 



£(fc,s) =r)(k;s,s). 



(103) 



Also in (|102|l . since both t and are arbitrary, such that 
we can make t ~ i', we may make the important assump- 
tion that 



V{k;s,t') = £(fc,s) 

= v( k ;s,s) = T)(k,s). 



(104) 



This tells us that in the case of the r](k;s,t') term, we 
need only concern ourselves with the on-diagonal terms 
|35| rj(k; s, s) — rj(k, s), which is a Markovian simplifica- 
tion. 

Going back to l|94|) . we can now write it as 
^D(k,j,\k-j\;t) 



= 1 - 



(vk 2 + vj 2 + v |k - j| 2 ) + v (k, t) + v (j, t) 



77(|k-j|,t) D(k,j,\k-j\;t), 



(105) 



which along with l|88|) can be used to evolve the memory 
time. 



A. Single-time Markovianized LET equations 

The final equations for the single-time evolution may 
now be summarised as: 



d_ 

dt 



2vk 2 ^j C(k,t)=2 J d 3 jL (k, j) x 



xD(k, j, |k - j| ; t)C (|k - j| , t) [C (j, t)-C (k, t)} 
= -2r){k,t)C (k,t) , (106) 

r,{k,t) = - fd 3 jL(k,j)D(k,j, |k-j|;i) x 

r<|k ' ' ilt) [C{j,t)-C{k,t)], (107) 



C(k,t) 



and 



d_ 

di 



D(k,j,\k-j\;t) 



= 1 - 



(vk 2 + vj 2 + v |k - j| 2 ) + V {k, t) + r}(j, t) 



77(|k-j|,t) D(k,j,\k-j\;t). 



(108) 



11 



These equations can be solved numerically with some 
suitable choice of initial conditions: 



C(k,t = Q) 



E(k,t = Q) 



(109) 



where E(k, t = 0) is an arbitrarily chosen initial energy 
spectrum, and 



D(fc,i,|k-j|;t = 0) = 0. 



(110) 



The last of these initial conditions follows from the defini- 
tion of D(k,j, k — j| ; t) in equation (87) and this in turn 
implies, from (|107f) . that r](k,t = 0) = 0, as is expected, 
because the cascade has not yet begun at t = 0. 

This set of equations is almost identical to those of 
the Test Field Model (TFM) : 34j], the exception being an 
extra term on the right hand side of (105) when compared 
with the corresponding TFM equation. As before, this 
extra term guarantees compatibility with K41. However, 
it remains to be seen how well the single-time LET theory 
performs when computed for the standard test problems. 



VI. CONCLUSION 

We have seen that our time-ordering procedures, as 
reported in pj have allowed us to tidy up some aspects 
of the LET theory. In particular, we have been able 
to derive a single-time form of the theory, which we 
have Markovianized so that it can be compared with 
well-known models such as TFM or EDQNM. Such a 
comparison will require numerical computation and this 
will be the subject of further work. However, we shall 
conclude here with some remarks about the role of the 
fluctuation-dissipation relation in the Eulrean two-time 
closures, DIA, SCF and LET theories. 

As we have seen, LET uses the FDR either to de- 
rive the response equation or to be used instead of a 



response equation. That is, with the second-order covari- 
ance equations for C(fc; t, t') and C(k; t, t) we can specify 
R(k; t, t') through the FDR and this gives us the requisite 
set of three equations. 

In contrast, SCF works with JEJ for C(k; t, t), the DIA 
response equation H20|) for R(k; t, t') and the FDR to cal- 
culate C(k; t, t'). Calculations based on these three equa- 
tions are known to agree quite closely with those for DIA, 
consisting of equations (|15fl . (|16fl and l|20|) . 

In the case of the DIA, one may test the idea of a FDR 
by introducing a modified response function R'(k;t,t'), 
such that 



R'(k;t,t') = 



Q(k;t,f) 
Q{k;t',t')' 



(111) 



This quantity plays no part in the calculation. However, 
at each stage, R can be calculated from the above re- 
lationship and compared with the actual DIA response 
function R at the same stage. It is this comparison that is 
the basis of the observation that the FDR is quite a good 
approximation at smaller wavenumbers but is less good 
in the dissipation range j^. However, such a comparison 
assumes that the DIA response equation is 'right' and the 
FDR is 'wrong'. In fact we know that the DIA response 
equation does not possess the correct behaviour at large 
Reynolds numbers and therefore cannot be a standard of 
comparison. In our view, the comparison of DIA with 
LET is a fairer test of the use of the FDR for turbulence. 

Lastly, our derivation of the FDR is correct to 
second-order in renormalized perturbation theory. Ac- 
cordingly it is an approximation, but no more an ap- 
proximation than the second-order covariance equations 
(|15fl and (|16|) . Therefore its use with these equations, as 
in the LET theory, is entirely consistent. Nevertheless, 
we should draw a distinction between this situation and 
that in microscopic equilibrium systems, where the linear 
form (|27|l holds to all orders in perturbation theory. 
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